Temperature and Magnetic Field Dependencies of Condon Domain Phase 
in Lifschitz-Kosevich-Shoenberg Approximation 
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The temperature and magnetic field behavior of non-uniform diamagnetic phase of strongly corre- 
lated electron gas at the conditions of dHvA effect is analyzed. It is shown, that in the framework of 
Lifschitz-Kosevich-Shoenberg approximation the magnetic induction splitting, as well as the range 
of existence of Condon domains, are characterized by strong dependencies on temperature, magnetic 
field and impurities of the sample. 

PACS numbers: 75.20.En, 75.60.Ch, 71.10.Ca, 71.70.Di, 71.25.-s; 71.25. He; 75.40-s; 75.40.Cx. 
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I. INTRODUCTION 



Oscillations of the thermodynamic quantities of elec- 
tron gas in magnetic field are the result of the oscillations 
of the density of states when successive Landau levels 
sweep through the Fermi level 0. At high magnetic field 
and low temperature the instability of strongly correlated 
electron gas due to magnetic interaction between elec- 
trons results in diamagnetic phase transition (DPT) Q 
with formation of Condon domain (CD) structure, which 
is extensively studied, both theoretically and experimen- 
tally 0-0. Magnetic interaction of strongly correlated 
electron gas results in non-linear dependence of local dia- 
magnetic moments on external magnetic field and tem- 
perature and gives rise such an exotic phenomenon as 
diamagnetic hysteresis 0], Q. 

Recent progress in experiments on observation of CD 
structure in Ag [f| and Be provides a natural stimulus 
towards a more detailed understanding of the properties 
of strongly correlated electron gas in the conditions of 
dHvA effect. 

Although the theoretical aspects of the CD forma- 
tion have been recently reviewed Q, some important 
questions concerning the DPT and non-uniform phase of 
strongly correlated electron gas are still open. In particu- 
lar, the detailed theoretical investigation of the influence 
of the temperature T, magnetic field (j,qH and impurities 
of the sample on such important characteristics of the 
non-uniform diamagnetic phase, as the range of existence 
of the CD structure and magnetic induction splitting, is 
still lacking. 

We present the theoretical investigation of the tem- 
perature and magnetic field behavior of the CD phase 
in the framework of Lifshitz-Kosevich-Shoenberg (LKS) 
approximation [lj , and compare our results with avail- 
able data on investigation of CDs in Ag 0],[H- 



II. MODEL 

The oscillator part of the thermodynamic potential 
density in the LKS approximation [2( can be written in 
reduced form with taking into account the shape sample 
effects: 



fi = 



47rfc 2 



a cos b + 2 a2 (l ~ n ) sm2 ^ 



(1) 



where b = k(B - (i Q H) = k[h ex + 4tt(1 - n)M] 
r • ; 1 n)y, /io-ff is the magnetic field inside the ma- 
terial k = 2ttF/(^oH) 2 , F is the fundamental oscillation 
frequency) , h ex — H ex — H is the small increment of the 
magnetic field H and the external magnetic field H exi 
n is the demagnetization factor. The validity of ther- 
modynamic potential density Q is restricted by applica- 
bility to homogeneous phase only, where the conception 
of demagnetization coefficient n is justified 0- In the 
conditions of strong magnetic interaction, when 



a(n H,T,T D ) > 1, 



(2) 
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a state of lower thermodynamic potential is achieved 
over part of dHvA oscillation cycle by the formation of 
CD structure. The usually observed diamagnetic domain 
structure is of stripe-domain type In the domain 

state, when the reduced amplitude of dHvA oscillations 
satisfies to the condition (|5J), the demagnetization factor 
n is replaced by the coefficient a 9j, which depends on 
magnetic field, temperature, impurities and geometry of 
the sample. 

The equation a(fioH,T,Tr)) — 1 defines the critical 
surface in three dimensions ^j,qH — T — Td . Above this 
surface the uniform diamagnetic phase exists, but below 
it, the CD phase appears in the part of the period of 
dHvA oscillations. 

In the non-uniform phase, when the condition is 
fulfilled, the metastable states, formed by magnetic field 



2 



h, are determined by 

an 

dx 



0. 



d 2 n 

dx 2 



>0, 



(3) 



identifying two local minima of the thermodynamic po- 
tential (JTJ). The diamagnetic moments y± corresponding 
to these minima, being the functions of magnetic field 
increment x and reduced amplitude of the dHvA oscil- 
lations a — a(/ioff, T, Td) 0, are characterized by the 
strong dependence on magnetic field, temperature and 
impurity of the sample and contribute to the magnetic 
induction splitting 5b = y + — ?/_ between two adjacent 
domains. The average magnetic induction splitting can 
be defined as the magnetic induction splitting at the cen- 
ter of the period of dHvA oscillations: 



Sb(x; a) « 6b(0; a) = y+(0; a) 



In this case the magnetic induction splitting 
lated due to 



Sb = 2a sin — . 



y-(0;a). (4) 
is calcu- 



(5) 



In limit a — ► 1 + + , e.g. in the nearest vicinity of the 
point of DPT, from the Eq. (JSJ) one can obtain the ex- 
pression for magnetization y± = ±y/[6(l — TJa)] 0- Ac- 
cording to Eq. JSJ magnetic induction splitting 8b = Sb(a) 
has strong dependencies on magnetic field fioH , temper- 
ature T and Dingle temperature Td through the reduced 
amplitude of dHvA oscillations a = a(/io-ff, T, To). 




The range of existence of the CD structure A in every 
period of dHvA oscillations can be estimated by consider- 
ation of the instability points, or inflection points, where 
the next two conditions must simultaneously hold 



on 



= 0, 



d 2 n 



dx dx 2 
Using the Eqs. (© one can obtain 



= 0. 



A = 2(y/~c 



1 



cos 1 — ). 

a 



(0) 



(7) 



From Eq. J7J it follows that the range of existence of the 
CD structure is also dependent on magnetic field (j,qH, 
temperature T and Dingle temperature Tjj. 

For proper calculation of the reduced amplitude of the 
dHvA oscillations a and investigation of the tempera- 
ture and magnetic field characteristics of CD phase, the 
correct topology of Fermi surface has to be taken into 
account. 

In case of ellipsoidal Fermi surface, the temperature 
and field dependence of the reduced amplitude of dHvA 
oscillations a is defined by 0,01 



a = a (^j, Q H) 



A(Ak)g,r) 
sinh X(^lqH, T) 



exp[-X(n H,T D )], (8) 



where 



2TT 2 kB'm c cT 



(9) 



m c is the cyclotron mass, fcs is the Boltzmann constant, 
e is the absolute value of the electron charge, c is the light 
velocity, h is the Planck constant, and Td = h/2irkBT is 
the Dingle temperature, which is inversely proportional 
to the scattering lifetime r of the conduction electrons. 
The limiting amplitude 



a = a(T — » 0,T £ 



0)=(^-) 3/2 (10) 



is the combination of temperature-independent factors in 
the Lifschitz-Kosevich formula , and 



FIG. 1: The magnetic induction splitting Sb as a function 
of the range of existence of non-uniform phase A. Closed 
circles show the different values of reduced amplitude of dHvA 
oscillations a: 1 - a = 1, 2 - a — 1.3, 3 - a = 2.6, 4 - a = 3.9, 
5 - a — 4.6, 2 - a — + oo. Both variables < 8b < 2w and 
< A < 2-7T are restricted by the period of the oscillations (2n 
in reduced units), which is shown by dash lines. The straight 
line between points 5 and 6 corresponds to slow increase of 
induction splitting Sb with increase of a in the interval a 6 
[4.6, oo), when A has reached its maximum value A = 2n. 
The open circle corresponds to the data Q . 



fi H m = (10.477ef0 2 / 3 (11) 

is the limiting field Q above which DPT does not oc- 
cur at any temperature T, tp is Fermi energy in eV, 
•q = m c /m and m electron mass. The validity of the 
LKS approximation and, consequently, the expression for 
reduced amplitude of the dHvA oscillations a © is re- 
stricted by the application to the spherical (or almost 
spherical) Fermi surface sheets, which is the case of noble 
metals |2| . The recent experimental data on observation 
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FIG. 2: Magnetic field dependence of the reduced range of 
existence of non-uniform phase A (a) and corresponding mea- 
sured range Acd (a) in silver for Td = 0.1K and four different 
temperatures: 1 - T = IK, 2 - T = 1.5K, 3 - T = 1.75K, 
4 - T — 2K. The straight horizontal line at (a) corresponds 
to the maximum value of A = 2-k. This line transforms into 
second power of magnetic field dependence of Acd = A/k 
at the same interval of magnetic field (b). At T = IK Acd 
reaches its maximum value at fioH « 15T. 



of CD structure in Ag |TT| , by measurement of the 
amplitude of the third harmonic of the ac susceptibility 
reveals the possibility to construct the phase diagram. 
As it was expected, the measured phase diagrams for Ag 
[lfij are in a good agreement with the calculated in the 
framework of LKS formalism phase diagrams, justifying 
the applicability of the Eq. (jHJ for belly oscillations in 
Ag. 



III. RESULTS AND DISCUSSIONS 



parameter a <E [1, 4.6] (Fig. [IJ. At a = 4.6 the reduced 
range of existence of CD structure spreads through all 
period of the dHvA oscillations and remains constant at 
further increase of a. 

The results of numerical calculations (Fig. show 
rapid increase of magnetic induction Sb with increase of 
reduced amplitude o in the vicinity of the point of DPT 
with slowing down of this increase at a S [2, 4.6]. Con- 
trary to it, the range of existence of the CD structure 
A is characterized by slow increase near the DPT point 
and rapid increase at the biggest values of a till its full 
value A = 2tt at a = 4.6. Thus, Sb = 2.43 >> A = 0.27 
at a — 1.3. At the interval a G [4.6, oo] the function 
Sb increases slowly with increase of a, while the range of 
existence of CDs remains constant A = 2tt. In famous ex- 
eriment on observation of CD structure in Ag by NMR 
the value of a was found to be about 2.6, which is con- 
sistent with the value of a, calculated from Eq. (jHJ at the 
conditions of the experiment: fi H = 9T, T = 1 AK and 
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When the condition @ is fulfilled, the sample is di- 
vided into domains with up and down magnetization. 
Both the magnetic induction splitting defined in explicit 
form by the Eq. © and the range of existence of CDs 
Eq. are the functions of reduced amplitude of the 
dHvA oscillations and cannot be considered as indepen- 
dent functions at wide range of applied magnetic field 
and temperature. Mathematically, the Eqs. |2Jl,© de- 
fine the function Sb — 56(A) in parametric form with the 



FIG. 3: Magnetic field dependence of the reduced range of 
existence of non-uniform phase A (a) and corresponding mea- 
sured range Acd (a) in silver for T = 1.2K and four different 
Dingle temperatures: 1 - To = 0.1K, 2 - Td = 0.25K, 3 
- T D = OAK, 4 - T D = 0.55JT . Similar to the Fig. H the 
straight horizontal line at (a) corresponds to the maximum 
value of A = 2n. This line transforms into second power of 
magnetic field dependence of Acd = A/k at the same inter- 
val of magnetic field in (b). Acd reaches its maximum value 
at noB w 20T. 
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Td = Q.8K. Using the reported values of the dHvA pe- 
riod AH « 17G, 5B/fj, ps 12G and evaluating the range 
of existence of the CD structure as Acd ~ 7G, one can 
estimate = 2ttSB/^ AH w 4.3 and A = kA C D ~ 2.6. 
These values are very close to 5& = 4.4 and A = 2.5, 
calculated from Eqs. (@J and JJJ). 

Due to the bell-like shape of the diamagnetic phase 
diagrams T = T(/j, H,T D ) (see, e. g. j8j) there is one 
critical temperature T c at given magnetic field and two 
critical values of magnetic field H T (H- < H + ) at given 
temperature. Another possibility for realization of the 
DPT is related to the concentration of impurities in the 
sample, which influence on the amplitude of the dHvA 
oscillations through the scattering r of conduction elec- 
tron. 

Near the point of DPT a -> 1 + 0+ (T -> T c + 0", or 
H — > iJzp + or Td — > Td. c + 0~) the temperature and 
magnetic field dependencies of magnetic induction 8b can 
be represented as follows: 



2[6A C T(A C )]V 2 (^_Z) 1 / 2 



, T^T c + 0" 



8b={ 2;<iAD, :■' -{±2__Jjly \ Tj „ . .. Tll , .+ .„■ 



Tr 



where L(t) = cotht — 1/t is Langevin function and 



(12) 



a T = ±[-- + A T L(A T )-A^]. (13) 

Here, A c = \(fi H,T c ), X D . C = \(/j,oH,T DiC ), A T = 
\(lioH Tl T), A^°' = \([j,qH t ,Td). It can be shown that 
near the point of DPT the range of existence of CD struc- 
ture A is related to the magnetic induction splitting 5b 
according to 



A - £<*)■. 



(14) 



It follows from Eqs. I|12l) and l|14l) that the parameters 
of the CD phase are characterized by the critical behav- 
ior near the DPT. The temperature and magnetic field 
dependencies of the order parameter 8b <|12fl can result 
in existence of soft mode. The possibility of softening of 
the orbital magnon mode was studied in 16]. 

In general case the temperature and magnetic field 
characteristics of CD structure can be calculated from 
Eqs. (@J and Q with taking into account the Eq. JSJ for 
reduced amplitude of the oscillations. To illustrate the 
behavior of the non-uniform diamagnetic phase, we cal- 
culated the temperature and magnetic field dependencies 
of the range of existence of CDs A and magnetic induc- 
tion splitting 8b for Ag, where the domains were observed 
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FIG. 4: (a) Temperature dependence of the range of exis- 
tence of CD phase for six different Dingle temperatures: 1 - 
T D = Q.IK, 2 - To = 0.25ff, 3 - To = OAK, 4 - To — 0.55K, 
5 - To = 0.7K and 6 - To = 0.85K. (b) Temperature 
dependence of the magnetic induction splitting due to CDs 
for four different Dingle temperatures: 1' - To = 0.1-R", 2' - 
To = 0.5K, 3' - To = 0.9K and 4' - To = 1.3JC. The strait 
horizontal line at (a) corresponds to the maximum value of 
A = 2tt. 
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and the DPT diagrams were confirmed experimentally 



The measured, or absolute, values of the range of exis- 
tence of CD structure Ac d and magnetic induction split- 
ting AB are defined as follows 



*CD 



A 
J 



(JM>H) 



2ttT 



2ttF 



(15) 
(16) 



The results of numerical calculation of the temperature 
and magnetic field dependencies of the parameters of the 
CD structure, A Eq. (JJJ and 8b Eq. (J5J, are illustrated 
in Figs. EE 

Fig. |21 shows the magnetic field dependencies of the 
range of existence of CD structure in Ag at constant Din- 
gle temperature Td = 0.1-fsT and several temperatures: 
T = IK, T = 1.5K, T = 1.75K, and T = 2K. The 
straight horizontal line in the graphical representation of 
the function A = A(^H) at T D = 0.1K (Fig. (a), 
curve 1) corresponds to the full value of reduced range 
A = 27r with the amplitude of oscillations a > 4.6 (com- 
pare with the straight line 5-6, Fig. 0. Corresponding 
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measured, or absolute, range of non-uniform phase exis- 
tence Ac d tEH is characterized by the second power field 
behavior at the same applied field interval (Fig.EJ. The 
function Acd has maximums at the middle range of ap- 
plied field [IqH 6 (15, 25T) and decreases with increase 
of the impurity of the sample. 

In Fig. |3| the range of existence on CDs is shown as 
a function of magnetic field for T = 1.2K and four dif- 
ferent Dingle temperatures: Td = Q.\K, Td = 0.25K, 
T D = OAK, and T = 0.55K. The growth of Dingle 
temperature results in decrease of the range of the non- 
uniform phase. Similar to Fig. [21 the straight horizontal 
line in Fig.[3Ja) corresponds to the maximum of reduced 
range A, which results in quadratic increase of the abso- 
lute range of CD existence Acd Q5t - 

The temperature dependencies of the range of exis- 
tence of CDs A (J7J) (Acd (|15|l 'l and magnetic induction 
splitting Sb © (SB JTSJl) at fi H = 10T and different 
Dingle temperatures Td are illustrated in Fig. 0] At the 
temperature interval T ~ — 1.5K for Td = 0.1K (curve 
1 (a)) the range of CD existence A spreads over all pe- 
riod of dHvA oscillations, reaching its maximum possible 
value. Near the point of diamagnetic phase transition, 
when a — > 1 + + the parameter of the order Sb exhibits 
critical behavior in accordance with Eq. (jl2|l . It results 
in softening of orbital mode. The last circumstance is im- 
portant for experimental observation of the DPT point 
and constructing the phase diagrams fioH — T — Td by 
measuring the temperature dependence of parameter of 
order AB = AB(T). The increase of the Dingle temper- 
ature Td results in decrease of the range 5 (Fig. Ufa)) 
and splitting Sb (Fig.gfb)). 

The magnetic field dependencies of induction splitting 
Sb (SB) are shown in Fig. © for five different tempera- 
tures and in Fig.[S]for four different Dingle temperatures. 
The maximum values of absolute induction splitting (SB) 
are shifted into high-field range [XqH ~ 25 — 40T. 

The magnetic induction splitting due to CD struc- 
ture changes the distribution of magnetic field in vac- 
uum above the surface of the sample. The possibility of 
detection of these changes by means of Hall probe spec- 
troscopy and mapping of CD structure is discussed in 
[12, O and realized in @, 0. It was found [H, [3 
for the plate-like sample in applied magnetic field along 
the normal to the surface H || Z with periodic domain 
structure along the Y-axis , that the normal component 
of magnetic field H n is described by the following peri- 
odic function of Y with the period 2D equal to the period 
of CD structure : 



Hn 



SB , cos(wY/D) 

tan" 1 —A — 

717x0 sinh (ttZ/D) 



(17) 



The measured value by method of Hall probes is the 
maximum splitting of the normal components of non- 
uniform magnetic field above two neighboring domains 
SH n , e. g. the maximum difference between two val- 
ues, H\ and H^, of normal component of magnetic field 
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FIG. 5: Magnetic field dependence of the reduced induction 
splitting Sb (a), caused by two adjacent domains, and corre- 
sponding measured value SB (b) in silver for Td = 0.1K and 
five different temperatures: 1 - T = IK, 2 - T = 1.5K, 3 - 
T = 1.75K, 4 - T = 2K, and 5 - T = 3K. The second power 
law dependence of the period of the dHvA oscillations on the 
applied magnetic field shifts the maximum values of measured 
characteristics to the high-field range. 



Eq. H17(l caused by magnetic induction splitting SB inside 
the sample 



SH n =Hl- Hi = 



2 SB 

7T fl 



■ tan 



sinh (ttZ/D) ' 



(18) 



The function l|18l) is of rapid decrease of Z. Thus, the 
magnetic field splitting above the sample Eq. Q16[l at the 
distance of one-third of the period of CD structure is 
about 0.155B / fiQ. It follows from Eq. I|18|) that the mea- 
surement of the magnetic field splitting SH n at given dis- 
tance Z = d above the surface of the sample can give in- 
formation about magnetic induction splitting SB inside 
the sample due to CD structure with given period 2D. 

The Hall probes measurements of mag netic field distri- 
bution above the sample surface [5j , |15| were performed 
on a high quality single crystal of Ag 2.4 x 1.6 x 1.0 mm 3 . 
The Dingle temperature was estimated to be Td = 0.2K. 
At the conditions of the experiment, [IqH = 10T and 
T = 1.3K ,5j, the magnetic induction splitting, cal- 
culated in LKS approximation with a Eq. (JSJ, gives 
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5B(1.3K) I [Iq = 16. 7G. These estimations are in agree- 
ment with the earlier famous NMR measurements in Ag 
[H , where at the slightly less favorable conditions for ob- 
servation of CD structure, e. g. fioH = 9T, T = 1AK 
and Td = 0.8K, the smaller value of induction splitting 
in the sample SB/^q k 12G was reported. At the condi- 
tions of the experiment [4| , the Eq. I|ltj|) with taking into 
account Eq. © gives the value of 5B(IAK)J fj, = 11. 7G 
close to the reported one, justifying the applicability of 
the LKS approximation 0] for belly oscillations in Ag. 

Correct treatment of the experimental results on mea- 
suring of distribution of magnetic field by Hall probe 
technic 5] can be done with knowledge of the ratio of 
the distance between the set of Hall probes and the sur- 
face of the sample d to the period of the domain structure 
2D (see, e.g. Eq. (|T%)l h Unfortunately, the distance be- 
tween Hall probes and sample surface was not reported 
in 0, and there are no accurate data about the char- 
acteristic magnetic sizes of the non-uniform phase, ex- 
cept of the statement, that the domain period was not 
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FIG. 6: Magnetic field dependence of the reduced induc- 
tion splitting 8b (a), caused by two adjacent domains, and 
corresponding measured value SB (b) in silver for T = 1.3K 
and four different Dingle temperatures: 1 - Td = 0.1-ft', 2 - 
T D = 0.5K, 3 - T D = 0.9K, and 4 - Td = 1.3K. Similar to 
Fig. 0, the second power law dependence of the period of 
the dHvA oscillations on the applied magnetic field shifts the 
maximum values of measured characteristics to the high-field 
range. 




T(K) 

FIG. 7: Comparison between theory and experiment, (a) 
Applied magnetic field dependence (a) and temperature de- 
pendence (b) of the magnetic field splitting 8H above the 
surface of the sample, caused by the induction splitting SB 
of CD structure. The solid lines correspond to the magnetic 
field distribution, calculated from Eq. 1161 in the framework 
of LKS approximation at the conditions of the experiment k|: 
HqH = 10T and Td = Q.2K. The circles are referred from [g. 
In calculation of the theoretical curves the fitting parameter 
d/2D = 0.1 is used. 



smaller than 150/im. The more complex measurements 
of the magnetic field distribution at different fixed po- 
sitions in the direction perpendicular to the surface of 
the sample would provide the lacking information about 
the period of the CD structure. To compare the theory 
with the available experimental results [5j we use the fit- 
ting parameter d/2D « 0.1, which gives reasonable value 
of d ks 15^tm for 2D = 150/im within accuracy of the 
experiment [l9j| . At d w 50/j.m which is the reasonable 
value for the average distance between unpolished surface 
and Hall probes, the calculated magnetic field splitting 
SH Eq. (|18|) is about 2.5G, which is at the edge of ac- 
curacy of the experiment |5( . It explains the absence of 
the signal from Hall probes for rough surface 5]. The 
results of comparison of the temperature and magnetic 
field dependencies of the measured magnetic field dis- 
tribution above the surface of the sample with the 
value of 6H n , calculated from Eq. I|18|) in LKS approxi- 
mation (JSJ) are present in Fig. 0. The reported critical 



7 



temperature T c w 3K at figH = 10T, defined from the 
temperature dependence of the measured splitting, is in 
accordance with the predicted value T c — 2.9K, calcu- 
lated in LKS approximation for Tjj = 0.2K at the same 
value of fi a H = 10T. 

IV. CONCLUSIONS 

We investigated theoretically the temperature and 
magnetic field characteristics of the CD structure in LKS 
approximation , • We show that the magnetic induc- 
tion splitting in non-uniform diamagnetic phase, caused 
by the arising instability of strongly correlated electron 
gas in high magnetic field and low temperature, and the 
range of the existence of this phase depend on temper- 
ature, magnetic field and impurity of the sample. We 
show that the magnetic induction splitting, e. g. the 
order parameter of the DPT, and the range of the exis- 



tence of this phase in every dHvA period are dependent 
on each other functions at the wide range of magnetic 
field with possible existence of the non-uniform diamag- 
netic phase. We calculated the critical behavior of the 
parameters of the non-uniform phase. The theoretical re- 
sults are in good agreement with available experimental 
data on observation of CDs in Ag. Further experiments 
on observation of electron instability and formation of 
CD structure would provide the lacking information on 
characteristic magnetic sizes, e. g. width of the domains. 
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